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Abstract 

The rich phenomena of the shear viscosity (77) to entropy density (s) ratio, n/s, in weakly 
coupled iV-component scalar held theories are studied, n/s can have a "double dip" behavior 
due to resonances and the phase transition. If an explicit goldstone mass term is added, then 
n/s can either decrease monotonically in temperature or, as seen in many other systems, reach 
a minimum at the phase transition. We also show how to go beyond the original variational 
approach to make the Boltzmann equation computation of n systematic. 



I. INTRODUCTION 



Scalar field theories are important tools to study spontaneous symmetry breaking. 
They can be used to demonstrate the Goldstone theorem by breaking a global symmetry 
and the Higgs mechanism by breaking a local symmetry. Furthermore, in some systems, 
scalar field theories are the low energy effective field theories of the underlying theories, 
such that model independent results can be obtained. A good example is the universality 
of the critical exponents for systems with the same symmetries near second order phase 
transitions. 

In the study of transport coefficients, scalar field theories also play an important role. 
It is proven that in a weakly coupled scalar field theory with quartic and cubic terms, 
summing the leading order (in the coupling constant expansion) diagrams in the Kubo for- 
mula for shear viscosity (77) is equivalent to solving the Boltzmann equation with effective 
temperature (T) dependent masses and scattering amplitudes 

Recently, there is a renewed interest in shear viscosity. It was conjectured {2] that no 
matter how strong the particle interaction is, 7]/s (jj per entropy density) has a universal 
minimum bound l/(4ir) in any system. This bound is motivated by the uncertainty 
principle and is found to be saturated for a large class of strongly interacting quantum 
field theories whose dual descriptions in string theory involve black holes in anti-de Sitter 
space [a-El. Much progress has been made in testing this bound and trying to identify the 
most perfect fluid with the smallest rj/s (see [t 9I] for recent reviews). It is found that r]/s 



can 



11 



3e as small as possible [10] (but still positive) in a carefully engineered meson system 



12], although the system is metastable. Also, in strongly interacting conformal field 



theories, 1/N r . corrections, with N c the size of the gauge group, can modify the rj/s bound 



slightly 



13|-ll6| 



In the real world, the smallest rj/s known so far be 



ongs to a system of hot and 



dense matter thought to be quark gluon plasma ( see |l7l-l20| for reviews) just above the 
phase transition t emp erature produced at RHIC 2lM29| with rj/s — 0.1 ± O.l(theory) ± 
0.08 (exp eriment ) 30]. A robust upper limit 77/s < 5 x 1/(Att) was extracted by another 



group 



31] and a lattice computation of gluon plasma yields rj/s = 0.134(33) [32]. Progress 



has been made in cold unitar y fe rmi gases as well. An analysis of the damping of collective 



oscillations gives rj/s > 0.5 33|, |34|. Even smaller values of rj/s are indicated by recent 



data on the expansion of rotating clouds 3Ja, [36] but more careful analyses are needed 



37]. 



In general, stronger interaction implies smaller rj/s and it is found that rj/s goes to 



a loca 



38 



minimum near the phase transition temperature (T c ) in a large class of systems 
-|40j. In particular, 7]/s develops a cusp(jump) at T c for a second(first) order phase 
transition and a smooth local minimum for a cross over. This behavior is seen in QCD 
with zero baryon chemical potential 



transition 



51 



38 



39 



4514501] and near 



;he nuclear liquid-gas phase 



521 ] . It is also seen in cold unitary fermi gases 



53|. in H oO, N, and He and 



in all the matters with data available in the NIST database [38|, [52|, [54|. Thus, it was 
speculated that this feature is universal. If this is indeed the case, then rj/s can be used 
to probe some parts of the systems which are hard to explore otherwise. For example, 
one can try to locate the critical point of QCD by measuring rj/s 52l. |55|. 

Theoretically, the behavior of rj/s going to a local minimum near the phase transition 
can be reproduced in controlled calculations of one-component weakly interacting real 



scalar field theories 401. In t 

n 

scalar field theories (see 141. 



lis paper, we extend the calculation to A-component real 



42] for earlier work in the symmetric phase and work in 



linear sigma models with QCD-like parameters [43|, l44j). Unlike the one component 
case, the resonance (called a) in the intermediate states makes the power counting more 
complicated. To perform a controlled calculation, we take the N — > oo limit and make 
the goldstone bosons (called 7r's) massive to stay away from the a pole. The resulting rj/s 
can decrease monotonically in temperature despite the phase transition. Also, when the 
resonance effect is important, rj/s can have a "double dip" behavior. It is conceivable that 
by tuning parameters, one can make the two local minimums to be close to each other and 
with the cusp at T c smoothed out such that one just sees the single minimum below T r as 



is shown in 



43j . These are behaviors different from that of the one-component theory 40 1 



and many other systems mentioned above. We also show how to go beyond the original 
variational approach to make the Boltzmann equation computation of r\ systematic. 



II. VALIDITY OF THE BOLTZMANN EQUATION 

n n 

It was proved first in weakly coupled </> theory |l| and then in hot QED [56j that the 
Boltzmann equation gives the sum of the leading order diagrams in the coupling constant 
expansion. Applying this equality to strongly interacting systems is dangerous because the 
Boltzmann equation might become invalid. As a semi-classical approach, the Boltzmann 



equation describes the evolution of the distribution function /(x, p, t), which specifies x 
and p at the same time. Thus, to be consistent with quantum mechanics, the size of the 

1 /2 

collision region (l sc ~ a sc , with o sc the scattering cross section) needs to be much smaller 
than the mean free path (l m f p ~ l/na sc , with n the number density). In other words, the 
quantum mechanical collisions happen in short distance (l sc ) while particles travel for long 
distances {l m f p ) as free particles between collisions. It is the wide separation of the two 
distance scales (7 sc // m /p *C 1) that allows the semi-classical treatment of quasi-particle 
collisions in the Boltzmann equation. 

For a theory with a dimensionless coupling constant g, if all the other mass scales are 
of the same order as the temperature T, then / sc // m / p ~ nal[ 2 oc g a , where a is positive 
for perturbation theories. For example, in a gas of massless pions, g ~ T/f n and a = 6. 
Thus the Boltzmann equation is applicable at low T Also, in one-component real 
scale field theories Q|, a = 3, with g the coupling of the four-scalar vertex. In strongly 
interacting systems, in general l sc /l m fp <C 1 is not satisfied and the Boltzmann equation 
is no longer applicable. 

In iV-component scalar field theories with a spontaneously broken O(N) symmetry, the 
u resonance could make l sc /lmfp ~ 1 even when the coupling is weak. We will first discuss 
the case with N = 0(1) then move to the large N case. In the symmetry breaking phase, 
the 7T7T scattering amplitude can be enhanced through the s-channel a pole in intermediate 
state. Adding the a — > ttti decay width T a which is O(X), the cross section a sc becomes 
O(X ) near the a resonance. This gives l sc /lmfp ~ 1 and the Boltzmann equation is no 
longer applicable. 

In the large N case, however, g can be scaled as 1/N such that the one- loop result is 
the same order in iV as the tree-level one. Therefore, l sc ~ g ~ 1/N, l m f p ~ l/(Ng 2 ) ~ N, 
and l sc /l m fp ~ 1/N 2 despite the a resonance. This is yet another example that large N 
systems are classical in nature. However, even though we will perform such a computation 
in the large iV limit, it is not clear whether l sc /lmfp "C 1 is sufficient to justify the use of 
the Boltzmann equation in the rj computation. 

One way to stay away from the a pole to obtain a reliable r\ is to add a mass term to 
7r's. Because if 2m n > m a , then a will never be on-shell in the tttt scattering. We will 
also present results in this case. 



III. MODIFIED O(N) MODEL IN THE LARGE N LIMIT 



We will study a iV-components scalar theory with the Lagrangian 

L = \(dj) 2 ~ \<tf - \m 2 ^ - ~ {$ 2 )\ (1) 

where = (7?, <Pn) and 7? = (4>i,(p2, • • • ,<f>N-i)- When m — 0, this theory has an 0(N) 
symmetry such that the Lagrangian is invariant under fa — > Rij<f)j with RijRji = 1. 
The m 2 term break the O(N) symmetry to 0(N — 1), thus it is called a modified O(N) 
model in this paper. As mentioned above, the inclusion of the m 2 term (m 2 > 0) is 
to avoid the production of on-shell a. a, b and m are renormalized quantities and the 
counterterm Lagrangian is not shown. The renormalization condition is that, at T = 0, the 
counterterms do not change the particle mass and the four-point couplings at threshold. 
The 1/N scaling in the b coupling makes sure that the quasi-particle masses are (9(-/V°). 
We will discuss the following cases: (I) a > 0, b > 0, the system is always in the symmetric 
phase. (II) a < 0, b > 0, in the m — > limit, the vacuum at T = breaks the O(N) 
symmetry spontaneously and there are N — 1 massless goldstone bosons. Those goldstone 
bosons become massive because of the m 2 term. At higher T, the symmetry is restored 
through a second-order phase transition. (Ill) Adding a term 



5L = — v NH4>n (2) 

to the Lagrangian of (II) to model a crossover. 

We will focus on the case of weak coupling in the large N limit and compute the 



effective potential via the standard Cornwall- Jackiw-Tomboulis (CJT) formalism [57] 
which has the one-particle irreducible diagrams included self-consistently. The detailed 
derivation of the finite N case is given in the Appendix A. We will only summarize the 
large N result here. 

In the symmetry breaking cases, we can expand shift the field = (7?, y/Nv + a) and 
expand the Lagrangian as 

£ = \ [(d^) 2 ~ <o* 2 + (d,af - ml fi a 2 ] - NU(v) 

—gi<j 3 — s^c?? 2 — Aicr 4 — A 2 <r 2 7f 2 — A 3 7? 2 7f 2 — A4O" , (3) 



where 



U(v) = -v 2 + h -v A + Hv, 
2 4 



rn 



l,o — a + m 2 + bv 2 , m 2 fi = a + 3bv 2 , 

bv b 
91 ~ 92 ~ ' 1 " 3 " 4iV ' 



A 4 = >VW • (4) 



Using the result of the Appendix A, the effective potential in the CJT formalism in the 
large iV limit reads: 

^P- = ^ 2 + + m +\f R ^ + ^ - 1] + \ Ll + 0(1/N) , (5) 
where L P = J R P(K,v) and P(Po) is the full (tree-level) propagators 

P~\K,v) = —K 2 + ml(v) , P -\K,v) = —K 2 + ml >0 (v) . (6) 
The condensate v = v is determined from minimizing the effective potential. It satisfies 

H = v [a + bv 2 + bL P ] + 0(1/ N) . (7) 
In case II, v ^ below T c , and 

m 2 n = — + m 2 + 0{l/N) , 

v 

m 2 a = ml fi (v) + bL P + 0(1/N) 

= — + 2bv 2 + 0(l/N). (8) 

v 

Note that, if if = and m — 0, the goldstone bosons (the n fields) remain massless below 
some critical temperature T c . The condensate v changes continuously to zero when T 
approaches T c from below showing that it is a second-order phase transition. 
Above T c , v — and 

m 2 = a + m 2 + bL P + 0(l/N) , 

= a + bL P + 0(l/iV) . (9) 

The entropy density of the system is given by the thermal dynamical relation s = 
-dV(v)/dT both below and above T c . 



IV. SHEAR VISCOSITY 



A. The Boltzmann Equation 

The equations needed to compute rj in the large N limit using the Boltzmann equation 
are derived in the Appendix B. The a distribution function is subleading, so only the 7r 
distribution functions are needed. All the 7Tj components are described by the same dis- 
tribution /""(x, p, t) = fp(x) (a function of space, time and momentum), whose evolution 
is described by the Boltzmann equation 



lr W*) = 7^7 / < 2 ™ {fif2 F 3 F p - F[F?fSf;} 



2N 



(10) 



12.-! 



where = 1 + ff, E p = a/p 2 + mf. The weighted measure i 



<^ri2;3p = |7i 



, (2%) 4 5 4 (k 1 + k 2 -k 3 -p)^j d 3 k t 



12;3p| 



n 



2'E 1 E 2 E 3 E r f \ (2ir)3 ' 
where T is the scattering amplitude for particles with momenta 1, 2 — > 3, p. 



11) 



ITi 



12;3p| 



26 + 



2„,2 



s — mi 



2b + 



Ab 2 v 



2„,2 



t — mi 



2b + 



4b 2 v 



2„,2 



u — mi 



(12) 



The first term is corresponding to the s-channel 7Tj7Tj — > TTjiTj scattering while the second 
and third terms corresponding to the t- and u-channel — > iTiiij scattering. As men- 
tioned in the introduction, the potentially divergent s-channel contribution is avoided by 
adding the 7r mass term such that the a would never be on-shell and the system remains 
a perturbative one. 

It is now straightforward to compute rj using the Boltzmann equation. We follow the 



same procedure as in 50|. It is known that computing rj in this approach is essentially a 
variational problem [53] . We will go one step further to show that the procedure we take 
can systematically approach the correct answer. Hence it does not rely on the Ansatz 
one takes in the computation. This procedure can be used in the computation of bulk 
viscosity as well. 



In local thermal equilibrium, the distribution function f p (x) = (e 



1)-' 



with f3(x) the inverse temperature and V^(x) the four velocity at the space-time point x. 
A small deviation of f p from local equilibrium is parametrized as 



f;{x) = f p (x) [1 - F» Xp (x 



(13) 



where F p =l + f. The energy momentum tensor is 



T ^= N Sw?it r * (x) - (14) 

We will choose the V(x) = frame for the point x. This implies d u V° = after tak- 
ing a derivative on V IJi {x)V^{x) = 1. Furthermore, the conservation law at equilibrium 
d/jT^lx^o = allows us to replace d t (3(x) and dtV(x) by terms proportional to V • V(x) 
and V/3(x). Thus, to the first order in a derivative expansion, Xpi x ) can be parametrized 

as 

Xp (x) = (3(x)A(p)V ■ V(x) + PWBMVMrtx) , (15) 

where i and j are spacial indexes, Bij(p) = B(p) (pipj — |<%) and VpV}] = 
iy iVj + V ' jVi — ■ V(a;)j /2. A and i? are functions of x and p, but we have sup- 

pressed the x dependence. 

Substituting f|T5|) into the Boltzmann equation, one obtains a linearized equation for 

B 

\PiPj ~ ^ijP 2 

= ^ J dT^T.TJliFl)-' [BM + B^ks) - B^h) - B ij (k 1 )\ , (16) 

where we have dropped the factor V[jVj] contracting both sides of the equation. There is 
another integral equation associated with V • V(x) which is related to the bulk viscosity 
C that will not be discussed in this paper. The V • and <9 t V terms in p^d^f will cancel 
each other by the energy momentum conservation in equilibrium mentioned above. 

In equilibrium the energy momentum tensor depends on pressure V(x) and energy 
density e(x) as T^(x) = {V(x) + e(x)} V^{x)V v {x) — V{x)5 lJ , u . A small deviation away 
from equilibrium gives additional contribution to T^ v whose spacial components define 
the shear and bulk viscosity 

5Tij = -T^V^x) + C«J y V ■ V(a;) . (17) 
5Tij can be computed using Eq.flT] 



The above two equations imply 



= {S\B). (19) 



Substituting Eq.fllEJ) into Eq.(JT9]) yields 

^ = 2uT / ^wZnVffiZKBM + 5 ^) - B^k 2 ) - Bij(kx)\ 

/ dT^FlFlTJl [Bij(p) + B tJ (k 3 ) - B l3 {k 2 ) - Z^)] 2 

</123» 



80T 



= (S|C|S), (20) 

where we have used that rj is invariant under 1 -H- 2, 3 -H- p, and 12 3p in the second line. 
It is easy to see from the more symmetric form in the second line that rj is non-negative 
so the matrix C in the third line is positive definite. 

B. A Variational Approach 

Now we review the arguments that the computation of rj can be formulated as a 



variational problem |45|, |58|. Let us rewrite Eq. f|T6|) as 

\S) = C\B), (21) 

whose projection onto \B) is just 

(S\B) = (B\C\B) (22) 

of Eqs. (|T9l and (|20l . Technically solving the projected equation (|22|) is easier than 
solving the integral equation f l2TT) . But this will give a wrong rj. However, 

77 = -(B\C\B) + 2(S\B) 

= -(B-C- 1 S\C\B-C- 1 S) + (S\C- 1 \S), (23) 

where \B - C~ X S) = \B) - C' 1 \S). Thus, if (|22J is satisfied but not (J2U, then r/ < 
(5 1 1 C 1 1 S) because C is positive definite. This implies that a variational calculation of rj 
is possible. One just demands ( )22|) and try to find the Ansatz that gives the maximum rj. 

C. Beyond Variation — Finding the Solution Systematically 

In a variational calculation, one starts with an Ansatz of B(p). Assume that B(p) is a 
smooth function, one can expand it using a specific set of orthogonal polynomials: 

f max 

B(p) = \p\yJ2brB ir \z(p)) . (24) 

r=0 



where z(p) = [3p, and B^ r \z) is a polynomial up to z r and b r is its coefficient. The 
overall factor |p| r will be chosen by trial and error to get the fastest convergence. The 
orthogonality condition 

1 



d 3 P | P | 2+ v 



f p FlB^\z)B^\z) = S^5 r , s 



(25) 



15T J (2vr) 3 E p 

can be used to construct the B^ r '(z) polynomials up to normalization constants. For 
simplicity, we will choose 

B<°>(z) = l. (26) 

Then, Eq. fpZO]) can be rewritten as 



Vtrial 



b) , 



(27) 



where \b) = (bo, b\ . . . , &- rmax ) and CV max is positive definite, while Eq. (fl9l can be rewritten 

as 

Vtrial = (S\b) . (28) 



Solving Eq.(j27P and Eq. (j25|l . we have |6> = C" 1 . 



S), and 



Vt 



rial — \ S C r ^ 



(29) 



Now, according to the orthogonality condition, rj = \S\b\ = NboS^°\ The other compo- 



nents of 



S) are zero, and 



oo 



(30) 



It can be shown that t] tria i increases with r max monotonically. Thus, one can approach 
the true r\ value systematically by increasing r max . The proof is as follows. C n+ i denotes 
a (n + 1) x [n + 1) matrix with elements = — n). Then the following identity 

holds: 



det (xy 



(31) 



'oo v " -/oo det (y ) det (Z) ' 

where X is a (n + 1) x (n + 1) matrix with elements Cy (i — 1 — (n + 1) , j — 2 — (n + 2)), 

Y is a n x n matrix with elements Cy = 1 — n), while Z is a (n + 1) x (n + 1) matrix 

with elements CV,- = 1 — (n + 1)). Both X and y are positive definite, so det (Y) > 

and det (Z) > 0. Furthermore, X is a real matrix, so det (X) 2 > 0. One then concludes 

that ( C~h ) < ( C~} 2 ) which implies that r] tria i increases with r max monotonically. 
V / oo V /oo 



Numerically, this algorithm converges very fast. For the case 11(a) shown in Fig. 1, 
using y = 1.89, r} tr iai increases by ~ 0.4%, ~ 0.03%, and ~ 0.003%, when r max increases 
from to 1, 1 to 2, and 2 to 3, respectively, for T = 30 — 300. 



D. Numerical Results 



In the large N limit, = (D(Nb~ 2 ) and rj = 0(iVV 2 ) by Eqs.flTBl). (!20l), and 
dT = 0(b 2 ). Combined with s = O(N), we have rj/s — 0(Nb~ 2 ). rj/ (sN) is shown in 
Fig. 1 for cases I-III. We see that the i]/s behavior could be different from that of the 
one component scalar model j^oj]. In case I, the system is always in the symmetric phase, 
and f]/s is monotonically decreasing in T. In case II, the system has a second order phase 
transition, rj/s decreases monotonically below T c , develops a cusp at T c . Then, depends 
on the parameters used, rj/s could be decreasing (11(a)) or increasing (11(b)) in T. In 
11(b), rj/s does not reach a local minimum at T c . We will discuss this case in more details 
later. Case III is similar to case II except that the cusp is smoothed out. 

In cases I-III, rj/s is monotonically decreasing near T = 0. This is because s approaches 
zero exponentially (7r's are massive) while rj approaches zero via power laws. This behavior 
persists even when 7r's are massless, but by a different reason. In this case, below T c , a can 
be integrated out. The resulting theory is a non-linear a model with massless goldstone 
bosons (the 7r's) that couple derivatively to each other. So 7r's become free particle at 
T = 0, and the interaction becomes stronger at higher T. Tow show this more explicitly, 
when m = H = 0, Eq.f lT2|) can be recast as 



Ab 2 



m 2 /s 



1 



Thus, 



m 2 /t 



T' 2 



1 



1 — m 2 j u 



r^Toc — 



(32) 



(33) 



with rricr approaching a constant at small T, the couplings between pions are weaker 
at lower T as mentioned above. And because smaller coupling implies larger rj/s, rj/s 
decreases monotonically near T = for massless tt's. 

As T — > oo, all the systems are in the symmetric phase. The only scales in the problem 
are m n and T. We find that rj/s has the 1/T expansion 

k N 



V 

s 



b 2 



ki a + m 2 

1 + 7t— + 



(34) 




FIG. 1: rj/(sN) vs. T for cases with a second-order phase transition (solid and dotted curves), 
a crossover (dashed curve), and with no phase transition (dash-dotted curve). The horizontal 
dotted line is the asymptotic line for rj/(sN) at high T, which is calculated from the massless 
case. Parameter can be in arbitrary units. 

with k ~ 571 + 80^ + 0(b) and fa ~ 0.84 - O.tfVb + 0(b). The leading term in the 
expansion is the straight line in Fig. 1 which corresponds to the i]/s for a theory with 
m = a = H = 0. (The slow running of the coupling b has been neglected. T is the 
only scale in the problem, so the dimensionless 77/s can only depend on the dimensionless 
coupling b but not T.) The leading T dependence comes from the a + m 2 term which has a 
positive sign for the symmetric phase (a > 0). For the symmetric breaking phase (a < 0), 
however, a + m 2 could still be positive. Numerically, this give a rj/s which does not reach 
a local minimum at T c for a second order phase transition as shown in case 11(b) of Fig. 
1. 

One might worry that whether case 11(b) is qualified as a second order phase transition. 
After all, the order parameter is defined on one component 0at whose mass is different 
from all the other N — 1 components. If we remove the 4>n, then the system does not 
have a phase transition in the first place. To answer this question, we study a similar 
model with just two real scalar fields js9|. One of the fields condenses below T c and the 
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FIG. 2: rj/(sN) vs. T for cases where the resonance effects are important (the thermal width 
of a is included). A "double dip" behavior is seen with the m = case. Parameter can be in 
arbitrary units. 

other stays in the symmetric phase. For simplicity, the interaction between the two fields 
is turned off. It is found that the rj/s behavior in this model is similar to that of 11(b). 



Our results of case II 



] a) and 11(b) differ from that of Ref. 43( which has the minimum 
of rj/s below T c . In [43||, the parameters of the Lagrangian are tuned to mimic the irn 
scattering in the real world. Thus, the S-channel tttt scattering diagrams are resumed 
to reflect the strong tctc scattering in intermediate energies. In our case, we keep the 
coupling b small and 2m 7r > m a to make sure irn scattering stays above the a resonance, 
such that we can apply the Boltzmann equation to compute rj reliably. If we set m = 0, 
such that m n = below T c , then the a can become on-shell in tttt scattering when 
T ~ 0(m a \T=o) — b x l 2 T c < T c . Then, indeed, a second local minimum (and sometimes 
also an absolute minimum, depending on the parameters) of rj/s below T c can be formed to 
have the "double dip" structure. Furthermore, it is conceivable that by tuning parameters, 
one can make the two local minimums to be close to each other and with the cusp at T c 



smoothed out such that one just sees the single minimum below T c as is shown in [43j. An 
rj/s plot with this feature is shown in Fig. 2 where the thermal width of o 60] is included 
in the computation of 77. However, the reader should be warned that while it might a 
generic feature to have a dip in rj/s below T c by a strong resonance, the 77 computed with 
Boltzmann equation in this case might not be reliable as discussed in Section 2. 



V. CONCLUSION 



We have discussed in details the computation procedure and the rich phenomena of the 
rj/s behavior in weakly coupled iV-component real scalar field theories. We have found 
that rj/s can have a "double dip" behavior due to resonances and the phase transition. 
It is conceivable that by tuning parameters, one can make the two local minimums to 
be close to each other and with the cusp at Tr smoothed out such that one just sees 
the single minimum below T c as is shown in [43j. If an explicit goldstone mass term is 
added, then rj/s can either decrease monotonically in temperature or, as seen in many 
other systems, reach a minimum at the phase transition. We have also shown how to go 
beyond the original variational approach to make the Boltzmann equation computation 
of rj systematic. 
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VII. APPENDIX A: CJT FORMALISM 



In this appendix we derive the effective potential of the Lagrangian 

£ = \{djf - if ■ $ - - \b (f) - l -c (f)' . (35) 

The if ■ $ term is included so one can mimic a cross-over with a non-zero H. The last 
term is a dimension six operator whose effect to rj/s is not studied in the main text but 



is included here for completeness. There could be two additional terms with dimension 
six: <p 2 ^0 • d 2 <pj, (j) ■ d 2 d 2 4> (the other terms are related to these ones via integration by 
parts). These terms can be removed by field redefinition or, equivalently, by applying 
the equation of motion. The inclusion of the dimension six terms shows that this is an 
effective field theory, which is valid under the cut-off scale 1/ (Ny/c) and is renormalized 
order by order in the momentum expansion Np^, p being a typical momentum scale in 
the problem, a, b, and c are renormalized quantities and the counterterm Lagrangian is 
not shown. The renormalization condition is that the counterterms do not change the 
particle mass and the four- and six-point couplings at threshold. 



We will use the standard Cornwall- Jackiw-Tomboulis (CJT) formalism [57] which has 
the one-particle irreducible diagrams included self-consistently. The effective potential in 
the CJT formalism reads: 
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where L s = J K S(K,v) and L P = J K P(K,v), and where S(S ) and P(Po) are the 
full (tree- level) propagators: 

S-\K,v) = -K 2 + ml(v),S^(K,v) = -K' 
P-\K,v) = -K 2 + ml(v) ,Pq X (K,v) = -K 2 + m 2 nfi (v) , (37) 

with the tree- level masses m 2 = a + 3bv 2 + 5cv 4 and m 2 = a + bv 2 + cv 4 . The expression 

I I 

for V is consistent with that of [61] in the c = limit. 

The self-consistent one- and two-point Green's functions satisfy 
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These yield 



H = v {a + bv 2 + cv 4 + (36 + 10cv 2 ) L s 

+ (N - 1) (b + 2cv 2 ) Lp + c (N 2 - 1) L 
+6c(N- l)L s L P + 15cL|} . 



(39) 
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(40) 



ml = ml Q + (b + Qcv 2 ) L s 

+ (N + 1) (b + 2cv 2 ) L P + c(N -1)(N + 3) L 2 P 

+2c (N + 1) L S L P + 3cL| . (41) 

In the large N, i.e. iV — )■ oo, limit, a sensible scaling is to make m a ^ = O(N ). This 
implies 

b = 0(N- 1 ) , c = 0{N~ 2 ) , 

v = O^N 1 ' 2 ) , H = O^N 1 ' 2 ) , (42) 

and V = 0(N). This is the scaling adopted in the main text. In this large N limit, 
the non-tadpole type loop diagrams are subleading in the effective potential. Thus, the 
Hartree approximation, which neglects the non-tadpole type loop diagrams, gives the 
correct result in the large N limit. For example, if H = 0, the goldstone bosons remain 
massless below T c as required. 



VIII. APPENDIX B: 1/JV EXPANSION OF THE COUPLED BOLTZMANN 
EQUATIONS 

In this appendix the leading contribution to the shear viscosity of an O(N) model in 
the large iV limit is derived. We start with the coupled Boltzmann equations: 
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and those for the other channels are defined analogously. 

Note that the 7Tj distribution is flavor independent. Thus, in Eq. ( I43p . there are only 
two independent distributions f 77 = f ni and f a . And the coupled Boltzmann equations 



can be written as: 
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where ^ = iV — 1 and g CT = 1. The scattering amplitudes (squared) are related to those 
in Eq. (JED as 
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In the large N limit, Eq.( l45|) is simplified to 



V£T\ 2 =\2 9 2\ 2 = 0(1/N) . (50) 



+jr [ dr^iftFZF^FZfU;} , (51) 

J 12 

^W) = +^ / 2 dr£T { av 2 -f; - , (52) 

where the iV dependence in <ir is factored out already, so all the N dependence is in the 
prefactors. The above equations imply 
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(53) 
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f / ^r^^(^)- 1 [^(p) - B£(A*) - BJj(ki)] , (54) 



5 
2 

where .By, defined in Eq.( IT5l) . describes how / changes when the velocity distribution is 
non-uniform. Eq.f l5"3"j) demands B^-{p) = 0(N). Eq. (f51|) demands B^(p) = O(N) such 
that Eq.([54D remains 0(1). Now, 



l) \B l ] = L f JP^fj?B l (p) (55) 

In the large N limit 

77 = NU [B w ] 

~ NL 77 [B n ] + cL a [B a ] , (56) 

where we have added a subleading term with prefactor c ~ O(N ). The final result for r\ 
should not depend on the choice of c. 

Substituting Eqs.( l53|54|) into Eq. (l56|) . one obtains 
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(57) 



Symmetries of the equations further gives 



n 



80 




(2 + c)/3 
60 



(58) 



By choosing c = —2, the subleading contribution can be subtracted. Thus, / and B^ 
decouple from 7] and the a contribution only appears in the intermediate states of tctt 
scattering. 

In summery, in the large N limit, one can use 



to solve Bf, and rj. 
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